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CN ! Abstract 

We determined the characteristic variability time scales (At p ) of 410 bright and 



long GRBs, by locating the peaks of their Power Density Spectra, defined and calcu- 
lated in the time domain. We found that the averaged variability time scale decreases 
CN ■ 

with the peak flux. This is consistent with the time-dilation effect expected for the 
^ ■ cosmological origin of GRBs. We also found that the occurrence distribution of the 

characteristic variability time scale shows bimodality, which might be interpreted as 
O ' that the long GRB sample is composed of two sub-classes with different variability 
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time scales. However, we found no difference for some other characteristics of these 
two sub-classes. 
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1. Introduction 

The isotropy of gamma-ray burst (GRB) directions and the deficiency of weak bursts 
suggest a cosmological origin of GRBs (Meegan et al. 1992), which has been confirmed by 
the detection of GRB afterglows (Costa et al. 1997; van Paradijs et al. 1997; Metzger et al. 
1997; Kulkarni et al. 1998). An energy release of ~ 10 51 — 10 53 erg (assuming no beaming), 
implied by the cosmological origin, makes them to be outstanding events in the universe. 

According to the most popular "fireball" model (Piran 1999), the observed 7-rays are 
emitted when an ultra-relativistic flow is converted to radiation. It has been suggested that 
the energy conversion occurs either due to an interaction with the external medium ( "external 
shocks") (Rees, Meszaros 1992) or due to collisions within the flow ("internal shocks") (Rees, 
Meszaros 1994). 

The progenitors of GRBs ( "central engine" ) remain the most mysterious aspects, since 
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they are hidden from direct observations. The most popular models are mergers of two compact 
objects (NS - NS merger or NS - BH merger), or the collapse of massive stars. Currently, it 
is difficult to distinguish them from observations. It is also likely that more than one type of 
progenitor could give rise to GRBs, and the overall GRBs may comprise different sub-classes 
corresponding to different progenitors. The most well-known classes are short and long classes 
separated at T90 ~ 2 s, arising from their bimodal duration (T90) distribution (Kouveliotou et 
al. 1993). In recent years, a few studies have been devoted to a reclassification of GRBs using 
multivariate analysis, which have yielded three sub-classes (Mukherjee et al. 1998; Balastegui 
et al. 2001). 

About 2000 GRBs 1 have been detected and recorded by the Burst And Transient Source 
Experiment (BATSE) on board the Compton Gamma- Ray Observatory. They have complicated 
and irregular time profiles which vary drastically from one burst to another. Recent studies 
concerning the internal shock paradigm show that the observed complicated temporal structures 
are directly associated with the activity of the central engine (Kobayashi et al. 1997). The 
observed variability provides an intersting clue as to the nature of GRBs (Piran 1999). 

Observationally, the light curve of a burst consists of successive pulses, which may be 
the fundamental emission units. Many people have decomposed light curves into pulses using 
model-dependent fitting (Norris et al. 1996) or peak finding selection (Nakar, Piran 2002; 
McBreen et al. 2001) and analysed the temporal properties of the resulting pulses, such as the 
rise-time, decay time, FWHM, and time interval between the neighbouring pulses. However, the 
real time profiles are so irregular that many pulses are overlapping and unseparable, prohibiting 
a clear decomposition of those pulses. 

In this work we adopted a different approach. We calculated the power density spectra 
(PDS) in the time domain for 410 bright long bursts, and regarded the time scale corresponding 
to the maximum of the PDS as the time scale of typical variations in the profile, and thus defined 
it as the characteristic variability time-scale of the burst. We then studied the distribution of 
the characteristic variability time scale and the correlation between the time scale with the 
GRB intensity. 

2. PDS in the Time Domain 

Many authors have calculated the PDSs of GRBs based on the Fourier transformation 
(Beloborodov et al. 1998, 2000; Pozanenko, Loznikov 2000). However, Fourier analysis cannot 
replace direct variability study in the time domain. Except for the periodic and quasi-periodic 
processes, there is no direct correspondence between a structure in the Fourier spectrum and 
the physical process taking place at a certain time scale. The power density at a given Fourier 
frequency can result from contributions from different processes on different time scales. 

1 ( http://gammaray.msfc.nasa.gov/batse/grb/catalog/current/). 
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Without using the Fourier transformation, a new technique for timing analysis in the 
time domain has recently been proposed (Li 2001; Li, Muraki 2002). Quantities characterizing 
the temporal properties, e.g., power density, coherence, and time lag, can be defined and 
calculated directly in the time domain with this technique. 

Following Li (2001), the variation power in a light curve, x(k), is defined as 

Var(x) ^k=A<k)-x] 2 
1 ' (At) 2 (At) 2 ' 1 ' 

where x(k),k = 1,...,N, is a counting series obtained from the time history of the observed 
photons with a time step of At; x = Y,k=i x (k) /N is the average counts. The power, P(At), 
represents the contribution from the variation at all time scales > At. The power density, p(At), 
can be defined in the time domain as the rate of change of -P(At) with respect to the time step, 
At. From two powers, P(Atx) and P(At 2 ), at two time scales, Ati and At 2 (At 2 > Ati), the 
power density at time scale At = (Ati + At 2 )/2 is evaluated approximately by 

For a noise series where x(k) follows the Poisson distribution, the noise power is given by 

_ Varjx) _ Jxl_ 

noiscl } ~ (At) 2 ~ (At) 2 ' 1 } 

and the noise power density at At = (Ati + At 2 )/2 is 

/a A _ -Pnoisc(Ati) — P no isc(At 2 ) , . 

Pnoisc[ } ~ At 2 -At x • [) 

The signal power density can then be defined as 

Psignal(At) = p(At) -p noise (At). (5) 

3. Data Reduction 

We used the BATSE Concatenated 64-ms Burst Data summed over energy channels II 
and III (50 — 300 keV), and we utilized the results of a background fitting for 64-ms Burst Data 
2 . In order to reduce the influence of the background, one would like to cut the light curve with 
a window neatly including the whole burst. We defined and used the "T 10 o window", where 
Tioo = Tgo/90%. The window was obtained by extending the Tgo window forward for 0.05Tgo 
and backward for 0.05T 90 , respectively. 

Since we are interested in using the PDS to investigate the variation power distribution 
over a large range of time scales, short bursts are not suitable for calculating the PDS, except 
that their durations are long compared with the timing resolution. We thus selected only bursts 
with T 10 o > 15 s. To avoid large Poisson fluctuations in the light curve, we excluded dim bursts 

2 ( http://cossc.gsfc.nasa.gov/batseburst/sixtyfour_ms/bckgnd_fits.html ). 
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Fig. 1. Time profiles of a burst (upper panel) and its PDS calculated in the time domain (lower panel). 
The solid line is the broken power-law model fit to the calculated PDS. 

with peak count rates < 250 counts per 64 ms bin. These two criteria gave a sample of 478 
bursts from the BATSE Current Catalog 3 . We then calculated the PDS of each burst in the 
time domain with the algorithm described in section 2. 

4. PDS Calculation 

For an individual burst, the power densities are calculated at time scales of At, arranged 
at equal intervals in logarithmic space. To avoid fluctuations in the PDS, we binned PDS with 
an approximately equal number of data points in each time-scale bin, and the uncertainties of 
the data were derived from the binning. 

Most of the bursts' PDSs showed a "bump" shape. The bump-shape indicates variation 
peaks at a specific time scale; we regard the peak time scale, At p , as being the characteristic 
variability time scale of the burst. 

To improve the accuracy of Ai p , we fitted the PDS with a broken-power-law model, 

p(At) = — — . - . — - — . - . - P 4 : the fit was acceptable for most bursts. The characteristic 

FV ; (At/Pi) p 2 + (At/Px) p 3 ' F 

variability time scale At p was determined from the peak location of the fitted PDS curve. 
Figure 1 gives an example of a burst's time profile and its PDS. 

3 ( http://gammaray.msfc.nasa.gov/batse/grb/catalog/current/ ). 
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5. Characteristic Variability Time Scales 



5.1. At p for a Simulated Time Series 

Does At p really reflect the variability time scale of the typical variation? How is At p cor- 
related with the sizes of pulses if we consider that a burst is composed of stochastic pulses? To 
answer these questions, we simulated a long stochastic time series with a 16 ms time resolution 
and 4800 s duration. The time series were divided into 100 segments, and the power densities 
at different time scales were calculated for each segment. The resulting PDS was obtained by 
averaging the PDSs over 100 segments. The pulse model was taken from Norris et al. (1996), 
with the following pulse shape: 

I(t)=Aexp[-(\t-t max \/r r )% t<t max , (6) 
= Aexp[-(\t-t max \/T d y), t>t max , (7) 

where t max is the time of the pulse's peak, A is the pulse amplitude, v is the "peakedness" 
parameter, and r r and Td are the rise and decay time constants, respectively. 

We assumed that the interval between the neighboring pulses is exponentially distributed 
with a mean value of 10 s, A is uniformly distributed within 10 — 200 counts (0.016s)" 1 , 
v = 1.2, the rise-time constants, 7~r ^ 3X6 cL Log-Normal distribution with a standard deviation 
cx(ln r r ) = 0.5, and the decay-to-rise ratio, r d /r r , is a constant. We defined the pulse rising 
time, t T , as the time during which the intensity, I(t), increases from 5%A to the peak. We then 
chose different values for the mean (In r r ) and the decay-to-rise ratio, t^/t t , to simulate the 
time series, and to calculate the PDS in the time domain. The characteristic variability time 
scale, Atp, was determined from the peak of the PDS. Figure 2 shows a simulated time series 
and the calculated PDS, respectively. Figure 3 plots At p vs. (t T ). Figure 3 shows that At p is 
truly the variability time scale of typical variations, and is sharply correlated with the rising 
time, t T , of the typical pulses. 

5.2. Distribution of At p 

Among the PDS samples, 63 PDSs keep rising untill the smaller time-scale limit where 
the PDS can be calculated, without showing any decrease. This means that the characteristic 
time scales of those bursts are smaller than, or equal to, the smaller time-scale limit, whereas 
they cannot yet be determined. After discarding the 63 samples and 5 badly fitted (x 2 / v > 50) 
samples, we obtained At p for 410 bursts. Figure 4 plots their histogram distribution; it is a 
bimodal distribution with a the demarcation at At p ~1 s. 

To investigate the reliability of the bimodality, we performed a chi-squared test (Press 
et al. 1992) for the differences between the distribution data and a known uniform distribution 
in the range of At p ~ 0.1 — 8.0 s, where 19 bins were included. The mean of the uniform 
distribution was set to that of the data sets. The test gave a chi-squared probability, Q(x 2 W) — 
1.437 x 10~ 4 , where the small value of Qix 2 ^) indicates a significant difference between these 
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Fig. 2. Simulated time scries (upper panel) and the calculated PDS (lower panel). The stochastic pulse 
parameters are (ln(r r )) = —2.0, Td/T r =2.5. 
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Fig. 3. Characteristic variability time scales At p vs. (i r ), for the simulated time scries. See definition of 
t r in the text. The solid line represents At p = (t T ). Note the correlation between At p and the average 
rising time, (i r ), of the stochastic pulses. 
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Fig. 4. Histogram distribution of the characteristic variability time scales for 410 bright long bursts. 

Table 1. Fitted results for time dilation. a 



Group 


Ai (s) 


Peak luminosity (10 50 erg s x ) 


■^min 


^max 


At p > 1 s 


2.30(±0.18) 


7.56(±2.81) 


0.28(±0.16) 


0.85(±0.32) 


At p < 1 s 


0.256(±0.004) 


5.52(±0.36) 


0.13(±0.11) 


0.72(±0.15) 



a The estimated uncertainties of Ato and the luminosity correspond to the 90% confidence level. 



two distributions. The bimodal distribution of At p indicates that bursts of At p > 1 s and At p < 
1 s may belong to two GRB sub-classes respectively. 

5.3. Time Dilation Test 

If a GRB occurs at a cosmological distance, then every structure in the time profile will 
be stretched by a factor 1 + z due to the expanding universe, where z is the red shift. Therefore, 
we should observe that the dimmer bursts have larger characteristic variability time scales than 
do the brighter bursts, assuming that bursts at different cosmological distances are "standard 
candles" with the same intrinsic characteristic variability time scale. 

According to the bimodal distribution of At p , we divided the bursts into the At p > 1 
s group and the At p < 1 s group; figure 5 is a plot of the distribution of the mean of At p 
in 5 brightness bins for each group, where the brightness is represented by P^, the peak flux 
measured at the 64 ms time scale. Both groups show that the averaged At p decreases with the 
brightness, which is consistent with the cosmological origin of GRBs. 

To quantitatively test the time-dilation effect, we fitted the At p vs. brightness distribu- 
tions with a model prediction; the model is described in the Appendix. Table 1 lists the fitting 
results for the At p > 1 s group and the At p < 1 s group, respectively, in a standard Q = 1, A = 
cosmology, with the normalized Hubble constant, h, set to 0.75. The luminosity is the product 
of the peak photon number luminosity and the assumed mean photon energy, 150 keV. The 
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Fig. 5. A< p vs. brightness distribution for the slow variable group (At p > 1 s) and the fast variable 
group (Atp < 1 s). Each bin includes an equal number of samples. The solid lines arc the best- fit model 
predictions. 

time dilation factors, (1 + z max )/(l + Zm in), derived for the At p > 1 s group and the At p < 1 s 
group are 1.46(±0.27) and 1.52(±0.20), respectively. 

6. Discussion 

We calculated the PDSs of 410 bright long bursts in the time domain, and determined 
their characteristic variability time scales by locating the peaks of their PDSs. The distribution 
of Atp is a bimodal distribution with the demarcation at At p ~ 1 s. GRBs may be divided 
naturally into two groups, a fast variable group and a slowly variable group. 

Following Kobayashi et al. (1997) and Piran (1999), in the internal shock paradigm, 
the emitted radiation from each collision between two relativistic shells will be observed as a 
single pulse, whose time scale depends on the cooling time, the hydrodynamic time, and the 
angular spreading time. In most of cases, the electron cooling time is much shorter than the 
latter two. The hydrodynamic time scale is determined by the time that the shock crosses the 
shell, whose width is d. A calculation revealed that this time scale (in the observer's rest frame) 
is of the order of the light crossing time of the shell, d/c. If the collisions of the shells take 
place at a larger radius, R, angular spreading (Sari, Piran 1997) affects the time scale of the 
pulse. If the distance between shells is 5, the resulting angular spreading time for the pulse is 
~ 5/c. If 6 > d, the observed variability time scale will be determined by the angular spreading. 
In any case, the observed temporal structure is directly associated with the variability of the 
"central engine", and the observed time scale is proportional to the size of the ejected wind 
or the separation between the winds. Consequently, it is likely that the observed time scale is 
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correlated with the size of the progenitor. 

Our results concerning the existence of a fast variable group and a slowly variable group 
may indicate that the bright long bursts comprise two sub-classes generated from different sizes 
of ejected winds or distances between winds, or from different sizes of the progenitors. Both 
groups show the trend that the characteristic variability time scale decreases with the brightness 
of the burst, consistent with the cosmological time-dilation effect. We fitted the data to a model 
for both the fast variable group and the slowly variable group, respectively. The smaller values 
of £ max and the time-dilation factors, compared with previous results (Norris et al. 1994; Che 
et al. 1997), may be because our sample does not include more weak bursts. 
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discussions and valuable advice. We thank S. N. Zhang for carefully reading the manuscript. 
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Appendix. Model Prediction for Time Dilation 

The cosmological time-dilation effect will stretch every time structure of GRBs if they 
are at cosmological distances. Assuming that GRBs have the same characteristic variability 
time scale, At , we should observe 

At = At (l + z). (Al) 

Since the red shift, z, is unknown, the observed peak flux, P [phs s -1 cm -2 ], is used to 
express z under the "standard candle" assumption. For simplicity, we introduce the peak photon 
number luminosity, L n [phs s^ 1 ], as the standard luminosity of GRBs, without considering the 
photon energy redshift. For a standard Friedmann cosmology (Q = 1,A = 0), the luminosity 
distance is 



d L = 2R (l + z- VTT^), (A2) 

where Rq = c/H = 9.25 h~ l x 10 27 cm, and h = H /100 is the normalized Hubble constant. 
Then, the observed peak flux is 

p = tV ( A3 ) 

Combining equations (A2) and (A3) yields 



0.3hJL n56 /P+ [l + 0.6hJL n56 /P 



1/2 



(A4) 



9 



where L n ^ = L n /10 56 is the normalized peak photon number luminosity. Therefore, from 
equation (Al) we obtain 



At 



1 + 0.3hJL n56 /P +11 + O.QhJL n56 /P 



1/2 



Air 



(A5) 



which can be compared with observations. 
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